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Abstract 



The influence of shape of expanding and rotating source on various char- 
pLn ' acteristics of the multifragmentation process is studied. The analysis is 

£\j ■ based on the extension of the statistical microcanonical multifragmen- 

CN , tation model. The comparison with the data is done for central Xe + Sn 

collisions at 50A ■ MeV as measured by INDRA Collaboration. 
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The multifragmentation process has been studied in a broad range of bom- 
barding energies and for various types of projectiles. The reaction mechanism 
is often considered in terms of two-step scenario where the first, dynamical 
step results in the formation of thermalized source which then, in the second 
step, decays statistically into light particles and intermediate-mass fragments 
(IMF's). Assuming that the thermal equilibrium is attained, various statistical 
multifragmentation models were employed for the second step (see [p], |[ and 
references quoted therein). These models were so successful that deviations 
between their predictions and the experimental data have been often taken 
as an indication for dynamical effects in the multifragmentation. However, 
one should be aware of several oversimplifying assumptions in the statistical 
calculations, such as , e.g., the spherical shape of the thermalized source. In- 
deed, one expects that the spherical shape is perturbed during the dynamical 
phase and the density evolution can give rise to complicated source forms [[3]. 
Even more important are the angular momentum induced shape instabilities 
M which may cause large fluctuations of both the Coulomb barrier and the 
surface energy even for moderately high angular momenta (L ~ AOh). In this 
paper, the non-spherical fragmenting source is considered within the statisti- 
cal model which is based on the MMMC method of the Berlin group [|IJ. 
The observables sensitive to the source shape are discussed and preliminary 
comparison with the experimental data for Eki n {Z) is presented. 

An explicit treatment of the fragment positions in the occupied spatial 
volume allows for a direct extension of the MMMC method to the case of 
non-spherical shapes ||. The source deformation in this case is considered 
as an external constraint, similarly as the freeze-out volume. Below we shall 
discuss axially symmetric ellipsoidal configurations (R x = R y ^ R z ) which are 
characterized by the ratio : 1Z = R x /R z . The freeze-out volume of deformed 
system is the same as that of an 'equivalent' spherical system with the radius 
R S ys = (RxRyRz) 1 ^ 3 - Consequently, the statistical weights in the Metropolis 
scheme |]I| due to the accessible volume for fragments remain the same. The 
Coulomb energy in our model is calculated exactly for every multifragment 
configuration of non-spherical nucleus. The general scheme to account for 
the angular momentum in the calculation of statistical weights is the same 
as discussed in Ref. @. The source deformation will change the moment of 
inertia of rotating system. In calculating the statistical decay of fragmenting 
system, the angular velocity of the source is added to the thermal velocity of 
each fragment. For each spatial configuration of fragments, part of the total 
energy goes into rotation and hence the temperature of the system will slightly 
fluctuate. We take also into account fluctuations of the moment of inertia 



arising from fluctuations in the positions of fragments and light particles. 

In calculating all accessible states in the standard MMMC method [|IJ, 
the source should be averaged with respect to spatial orientations of its axes. 
Some of these states will not be accessible if the angular momentum is con- 
served |6|, [7j . We disentangle in the MMMC code the beam direction (the z - 
axis) and the rotation axis (the x - axis, perpendicular to the reaction plane). 
The rotation energy is then : L 2 /2J X = L x 2 /2J X , where J x is the rigid-body 
moment of inertia with respect to the x axis. Averaging over polar angle 9 
is not consistent with the angular momentum conservation. On the contrary, 
averaging over 2ir in the angle corresponds to averaging over azimuthal an- 
gle of the reaction impact parameter and should be included. Averaging over 
rotation angle ip around L depends on the considered relationship between a 
rotation time : r rot = J x /L x , and a characteristic life-time of the source r c . 
For high angular momenta when r rot <C t c , the full averaging in < ip < 2tc 
should be performed. In the opposite limit when r rot 3> r c , only states with 
x/j ~ are accessible. 

In central HI collisions, part of the excitation energy can be stored in the 
compression energy of pre-formed source which is transformed into the kinetic 
energy of fragments during the collective expansion . To get some insight into 
the influence of collective flow on the multifragmentation process, we shall 
mimic this effect by adding the blast velocity v^ to the thermal velocity of 
each particle/fragment for any event simulated by the Metropolis method. 
The average collective energy of expansion, similarly as the deformation en- 
ergy, is not included in the value of the total excitation energy. We assume 
a simple scaling solution of non-relativistic hydrodynamics which yields the 
radial velocity profile [RL El : 



v b (r) = v (r/R o r , (1) 

where v and R are the strength and scale parameters respectively, and the 
exponent a characterizes the power-low profile function. R corresponds to the 
system size at the beginning of the scaling expansion regime. Hence, Rq < R sys 
and we take R = 0.7R sys in all studies. Strictly speaking, the profile function 
(HD describes radial expansion of a spherical source. For axially symmetric 
expansion, the velocity profile may be more complicated. But even in this 
case, the scaling solution ([[]) with 0.5 < a < 2 was successfully applied for 
describing the velocity profile of the transverse expansion in high energy heavy- 
ion collisions |9|, pi], |ll|]. In the multifragmentation case, we are dealing with 
a non-spherical expansion of unstable matter and therefore a may be treated 
as a free parameter. 



Let us now consider the multifragmentation of l97 Au having the angular 
momentum L = 40h and the thermal excitation energy 6 A ■ MeV. These 
parameters characterize an equilibrized source formed in central Xe + Sn col- 
lisions at 50 A- MeV studied by INDRA Collaboration fl2j|13l. All calculations 
are carried out at the break-up density p fa po/6, what gives R sys = 12.8 fm 
for the effective radius of 197 Au source. In the following, we consider ellipsoidal 
forms with the ratio of axes 1Z = 0.6 (the prolate shape) and 1Z = 1/0.6 = 1.667 
(the oblate shape). We have found that none of the observables related to the 
fragment-size distribution is sensitive to the source deformation at these high 
excitation energies. On the contrary, the c.m.s. angular distribution of the 
largest fragment [Z = Z max ) is a highly sensitive observable (see Fig.l). In 
the absence of collective expansion, the angular distribution of Z max is isotropic 
for oblate configurations and has small forward - backward peaks for prolate 
configurations if the ip- averaging is performed in the whole available interval : 
< ip < 2tt. For the 'frozen' spatial configuration (ip = 0), the shape effect is 
clearly seen: for the prolate form one finds strong forward - backward peaks, 

while for the oblate form the heaviest fragment is predominantly emitted 
in the sideward direction (8 cm = 7r/2), like in the hydrodynamic splashing 
scenario. The collective expansion (a = 2) enhances the 'deformation effect'. 
One may notice a strong increase of forward and backward peaks in the prolate 
case and the appearance of a strong peak at 9 cm = ti/2 in the oblate case. 
Similar features can be seen also in the cumulative angular distributions of 
all IMF's but the relative amplitude of the 'deformation effect' in that case is 
smaller. 

Large sensitivity to the source shape can be seen in the analysis on an 
event-by-event basis using global variables, such as the sphericity, coplanarity, 
aplanarity and the flow angle Qfi ow ||. In the latter case, the shape differences 
manifest themselves even for v^ — 0. The whole effect is extremely sensitive 
to the collective expansion and is enhanced furthermore for the 'frozen' con- 
figuration (ip — 0). It should be mentioned that ©/^-distribution is the only 
observable which explicitly depends on the angular momentum ||. 

We have compared the model predictions with the experimental data for 
central Xe + Sn collisions at 50 A ■ MeV [p], [E| [L| f0|. Various statistical 
multifragmentation codes assuming the spherical source shape (the standard 
MMMC 0, the SMM 0) have been tried earlier to explain these data. All 
of them were successful in predicting fragment partitions 0, [M], |j| which, 
as stated above, do not provide a constraint on the source shape. On the 
other hand, these models failed to reproduce the kinematical observables such 
as the IMF's angular distributions, the event shapes (sphericity, coplanarity, 



aplanarity), the IMF's average kinetic energies Eun(Z). In the latter case, 
the behaviour for heaviest fragments could not be reproduced even including 



radial expansion [12]. To understand this issue, we have applied our model 



for different shapes, orientations, expansion profiles (following (1)), starting 
from the source of 197 Au having 6 A ■ MeV of thermal excitation energy and 
an angular momentum of L = 40h. All calculated events have been filtered 
with the INDRA software replica, and then selected with the experimental 
centrality condition : complete events (i.e., more than 80% of the total charge 
and momentum is detected) and Qfi ow > 7r/3. We consider first the average 
kinetic energies of IMF's as a function of their charge. Results for different 
deformations 1Z = 0.6, 1.667 and expansion parameterizations are compared 
in Fig. 2 with the experimental data [12, O, H[. The calculations have 



been done for the 'frozen' configuration (ip=0). One can see that the data 
excludes the radial velocity profile with a=l and clearly favours the expanding 
prolate source with ct=2 and ip ~0. These source characteristics have been 
confirmed also by the observables related to the event shape, like the sphericity, 
coplanarity and aplanarity distributions, for which an excellent fit of the data 
has been obtained. Finally, we have examined the experimental IMF's and 
heaviest fragment angular distributions, as well as the ©/^-distribution which 
again clearly discriminate between prolate, spherical and oblate source shapes. 
On the basis of this whole set of observables, we believe to have the strong 
evidence in central Xe + Sn collisions at 50 A ■ MeV for the formation of 
prolate, slowly rotating source with the orientation of the source aligned with 
the beam axis. 

In conclusion, we have made an extension of the Berlin statistical mul- 
tifragmentation code in order to include the effects of deformation, angular 
momentum and collective expansion. Due to the change in the Coulomb en- 
ergy for deformed freeze-out configuration, the 'deformation effect' is clearly 
seen in the IMF's angular distributions (Z max - angular distribution) as well 
as in the Qfi ow - distribution. A surprising interplay between effects of non- 
spherical freeze-out shapes and the memory effects of nonequilibrium phase of 
the reaction, such as the rotation and the collective expansion of the source, 
has been found. The influence of shape on rotational properties of the system 
is not only reduced to the modification of the momenta of inertia. The limits 
on the averaging interval over the angle i/j about the rotation axis, which are 
defined by the time scales involved, affect strongly the angular observables 
and are able to enhance strongly the 'deformation effect'. These constraints 
may be important for certain observables used in experimental procedures of 
selecting a specific class of events. Other striking finding is that the collective 



expansion allows to disclosure the source shape in the analysis using global 
variables as well as in the study of Z- dependence of E^ in . All these experi- 
mental observables, including Eki n (Z), are well reproduced assuming strongly 
deformed (1Z = 0.6), slowly rotating (L = AOh) and expanding source which 
has the radial velocity profile (1) with a ~ 2. Alongside with the observables 
discussed here, it would be interesting to investigate the velocity and angular 
correlations between fragments which are sensitive to the source shape at the 



freeze-out. Such a work is now in progress [16 



Acknowledgements 

We are grateful to D.H.E. Gross for his encouragement and interest in this 
project. We also thank G. Auger, A. Chbihi and J. -P. Wieleczko for their 
interest in this work. 

References 

[1] D.H.E. Gross, Rep. Prog. Phys. 53 (1990) 605. 

[2] J. P. Bondorf et al, Phys. Rep. 257 (1996) 133. 

[3] K.K. Gudima et al., J. Phys. G: Nucl. and Part. Phys. 5 (1976) 229; 
W. Bauer et al., Phys. Rev. Lett. 69 (1992) 1888; L.G. Moretto et al., 
Phys. Rev. Lett. 69 (1992) 1884; B. Borderie et al., Phys. Lett. B 302 
(1993) 15. 

[4] S. Cohen and W.J. Swiatecki, Annals of Physics (N.Y.) 82 (1974) 557; 
M.E. Faber et al., Acta Phys. Pol. B15 (1984) 949. 

[5] A. Le Fevre, M. Ploszajczak and V.D. Toneev, |nucl-th/990i09£ . 

[6] D.H.E. Gross, Phys. Rep. 279 (1997) 119. 

[7] A.S. Botvina and D.H.E. Gross, Nucl. Phys. A 592 (1995) 257. 

[8] J.P Bondorf, S.I.A. Garpman and J. Zimanyi, Nucl. Phys. A 296 (1978) 
320. 

[9] T. Csorgo, B. Lorstad and J. Zimanyi, Phys. Lett. B 338 (1994) 134. 

[10] A. Polleri, J.P. Bondorf and I.N. Mishustin, Phys. Lett. B 419 (1998) 19. 

[11] K.S. Lee and U. Heinz, Z. Phys. C43 (1989) 425; Z. Phys. C48 (1990) 
525. 

[12] R. Bougault et al. (INDRA Collaboration), Proc. XXXV Int. Winter 
Meeting on Nuclear Physics, ed. by I. Iori , Bormio, (1997), p. 251. 



[13] N. Marie et al. (INDRA Collaboration), Phys. Lett. B 391 (1997) 15. 
[14] A. Le Fevre, PhD-Thesis, GANIL T 97 03. 
[15] S. Salou, PhD-Thesis, GANIL T 97 06. 
[16] A. Le Fevre et al., in preparation. 



1.5 






n 0.5 




o 




1.5 



0.5-- 



(a) 



(c) 




i i i i i i 



1 -0.5 0.5 1 

cos(0 



(b) 



(d) 



0.5 0.5 1 
) 



cm Zmax 

Fig. 1. Angular distribution of Z max in the c.m.s. Plots on the l.h.s. correspond to 
averaging over the whole available interval < ip < 2ir, whereas plots on the r.h.s. 
correspond to the 'frozen' configuration ip = 0. (a),(b) : L = 40/i, vu = 0; (c),(d) : 
L = 40^, v b = 0.08c, a=2 . 
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Fig. 2. Kinetic energy of IMF's is plotted as a function of Z. Thick lines correspond 
to events satisfying experimental 'centrality' condition : Qfi ow > ir/3 and accepted 
by the INDRA filter. Thin lines correspond to unfiltered events and no condition 
on Qfiow Prolate freezout configuration is denoted with full lines and oblate con- 
figuration is denoted by the dashed lines. The experimental data [13] are shown by 
dots, (a) L = 40h, a=2, vu = 0.08c. (b) The same as in (a) but for q=1. (c) Com- 
parison between the experimental data and the calculation (a=2) for both prolate 
and oblate freeze-out configurations and the blast velocities : v^ = 0.088c (prolate 
shape), Vfj = 0.063c (oblate shape), chosen as to reproduce the experimental mean 
kinetic energy per nucleon of IMF's : 5.2 ± O.lMeV/nucl. (d) The same as in (c) 



but for a=l. The blast velocities are : v^ - 
oblate freeze-out configuration respectively. 



0.08c and v^ = 0.07c for prolate and 



